Parametric Luminescence of Microcavity Polaritons 
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The spectral and dispersive emission properties are analytically determined for the two- 
dimensional system of exciton-polaritons in microcavities excited by a resonant and coherent optical 
pump. New collective excitations result from the anomalous coupling between one generic polariton 
state and its idler, created by the scattering of two pumped polaritons. The corresponding para- 
metric correlation is stimulated by the emitter and idler populations and drives very efficiently the 
luminescence. The intrinsic properties of the collective excitations determine a peculiar emission 
pattern. 



In their pioneering experiments, Weisbuch et al. |l]] 
discovered a new kind of two-dimensional quasi-particlcs, 
resulting from the strong coupling between quantum well 
excitons and confined photons in a semiconductor micro- 
cavity. These peculiar particles, called microcavity po- 
laritons, have a very sharp dispersion due to the very 
light mass of the cavity photon. Remarkably, this rep- 
resents a condensed matter system of small mass quasi- 
particles which can be manipulated through laser beams 
both in frequency and momentum space. In principle, 
the low polariton density of states could allow large oc- 
cupation numbers at relatively small densities of parti- 
cles , well below the critical saturation value due to the 
fermionic nonlinearities. In other words, the polariton 
system could exhibit bosonic properties || . Furthermore, 
unlike unbound electron-hole pairs in ordinary semicon- 
ductor lasers, microcavity polaritons have a relatively 
short time of recombination into external photons. All 
these ingredients are indeed very promising for applica- 
tions in the domain of ultrafast all-optical switching and 
amplification. 

Recently, Dang et al. || measured photoluminescence 
spectra from a II- VI microcavity excited with a nonres- 
onant pump. A threshold was observed in the depen- 
dence of the polariton luminescence intensity as a func- 
tion of the input power. Similar results were reported 
by Senellart and Bloch in a III-V microcavity [Q . These 
experiments have been interpreted in terms of enhanced 
scattering of reservoir excitons into the emitting polari- 
ton modes. The origin of the enhancement has been at- 
tributed to bosonic stimulation due to final state occu- 
pation. 

More recently, great insight into the subject has been 
given by angle-resolved experiments under resonant ex- 
citation l^-g. In this kind of experiments, polaritons 
are optically excited at a desired energy and momentum, 
allowing a direct control of the polariton dynamics. In 
particular, Savvidis et al. ^ have uncovered a new kind 
of polariton parametric amplifier through pump-probe 
experiments. The angular selection rules for the para- 



metric conversion of pumped polaritons into the probe 
and idler modes are unambiguously given by the energy 
and momentum conservation for the polariton scatter- 
ing. In the case of an applied probe, the nonlinearity can 
be explained in terms of phase-matched wave-mixing of 
polariton matter beams ||. Remarkably, giant nonlin- 
earities occur also when the probe beam is absent and 
the spontaneous emission is detected ||,|,|. Spectrally, 
the luminescence shows a very peculiar distortion with 
respect to the bare polariton dispersion HI). Moreover, 
the emission is characterized by a finite angular width 
which appears to be an intrinsic property of the system 
HU ■ Indeed, the whole phenomenology suggests that 
the interacting polariton matter provides a new kind of 
collective excitations. The nature of these excitations 
and the connection with the stimulated scattering mech- 
anism are not yet known. The explicit solution of this 
intriguing problem is the motivation and main result of 
our present investigation. 

In this paper, we present a theoretical analysis of the 
polariton nonlinear emission in the case of resonant and 
coherent pumping. The coherence induced by the pump 
field allows an anomalous coupling between a generic po- 
lariton state and its idler, originated by the fission of two 
pumped polaritons. The anomalous coupling gives rise to 
new collective excitations. The luminescence is driven by 
the emitter-idler parametric correlation which is stimu- 
lated by the emitter and idler populations. We establish 
how the intrinsic properties of the collective excitations 
determine the rich emission features. In particular, (i) 
the peculiar dispersion, (ii) the lineshape features and 
(iii) the two-dimensional pattern are provided by our re- 
sults. 

In order to investigate the nonlinear emission of mi- 
crocavities in the strong exciton-photon coupling regime, 
we work directly in the polariton basis. The polari- 
ton states have a two-dimensional character and are de- 
scribed by the in-plane wave-vector and their spin. In 
the following, each k will be a two-dimensional vector 
in the quantum well plane. Moreover, we will consider 
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the case of a circularly polarized pump beam. Neglect- 
ing the spin relaxation, only the polariton states with 
a definite circular polarization will be retained. When 
showing numerical results, we will use the angular co- 
ordinates (0 x ,6y) defined by (k x ,k y ) — ~(sm6 x ,sin9 y ), 
where u> is the emission frequency. Actually, these co- 
ordinates represent the angles of the far-field images in 
current experiments f§,§U§. Since we consider the situa- 
tion of resonant pumping of the lower polariton branch, 
we can neglect the nonlinear contribution due to the up- 
per branch. The destruction operator for a lower polari- 
ton with in-plane wave- vector k is p^ = Xkbk + CfeOk, 
where &k and ak are the exciton and photon Bose opera- 
tors respectively and Xk, Cfc the corresponding Hopfield 
coefficients (Xk > and Ck < 0). The lower polariton 
Hamiltonian g is H = H LP + Hp 1 / + H qm . The free 
term H^p = Ei,p(k) p\j)k contains the lower polari- 
ton energy dispersion Ei,p(k). The polariton-polariton 
interaction term reads 

jjeff _i V 3lv pp J J n, n, , 

PP ~~ 2 2^1 A V k,k',q Pk+q iV-q Pk Pk', 
k,k',q 

where Ax is the two-dimensional exciton radius and A is 
the macroscopic quantization area. The interaction po- 
tential is V££ A = {2^-(|C k+q | X v + |C k ,|X k+q ) + 

f|Ak + qAk4Ak'-qAk, with n sat = 7/(16ttA x ) the ex- 
citon saturation density and e the quantum well dielec- 
tric constant. Notice that Vjf^ is positive and rep- 
resents a repulsive interaction. The cavity system in- 
teracts with the external electromagnetic field through 
the standard quasi-mode coupling Hamiltonian H qm = 

J d£l{J2 k gCk p\. «k,o + H.c.} . The operator a k ,o ob- 
structs an external photon with in-plane wave-vector k 
and frequency tt. The spectrum of the polariton lumi- 
nescence is 

r+oo 

PL(k,t,u) oc \C k \ 2 TZe / dre^"-^ (p u (i+r)p k (t)) , 
Jo 

where Pk(t) is the time-dependent polariton operator. 
Here and in the following we use the Heisenberg repre- 
sentation of time-dependent operators and take expec- 
tation values on the stationary state. In the steady- 
state regime, the two-time correlation (p k (i + 7~)p k (t)) 
depends only on the relative time r, and not on t. 
This implies that the stationary spectrum can be cal- 
culated at a fixed time t = to when the steady-state 
regime is achieved. For simplicity, we choose to = 
and hence PL(k,ui) cx |Cfc| 2 7^e{(p k (w)p k (0))}, where 
PkH=/ +O °dre^ 0+ )>k(T). 

An applied cw-optical pump with in-plane wave- 
vector k p drives a polariton polarization (p k (i)) — 
e~ tUpt \(pk )\, where ui p is the laser frequency reso- 
nant with the lower polariton mode. We are in- 
terested in the generic polariton operator Pk(t) with 



k =/= k p . The polariton-polariton interaction cou- 
ples Pk(t) to the pumped mode and to the idler po- 
lariton operator p kidler (i) , where kidier = 2k p — k. 
This corresponds to the fission process {k p ,k p } — > 
{k, kjdier}- Disregarding the pump noise, we can 
perform the replacement P k (0P kidler (*)Pk p (£)p kp (i) — 

Pk(*)Pk id i er (*)(? , k p (i)) 2 - Furthermore, we neglect multi- 
ple scattering, that is interaction between modes other 
than the pumped one (i.e. {k, k'} — ► {k + q, k' — q}). 
The limit of validity of this approximation will be dis- 
cussed later on. With these assumptions, the Heisenberg 
equations of motion read 

;t . d ( Pk(t) \ _ fjpar ( Pk(t) \ 



-K (*)e~ £ 

Kidler V ' 



where the coupling matrix is 




FIG. 1. Lower Polariton energy dispersion Elp — -Eexc(O) 
(meV) versus the angles 6 X and 6 y (deg). A sketch is shown 
for a parametric fission of two pumped polaritons. 



The off-diagonal term contains the energy E£ lt = 
(^C P k p ,k-k p + ^k i p P kp,k P -k)/2 and the pump-induced po- 
larization Vk p = -^j(Pk p )- We point out that |"Pk p | 2 = 
n p \ 2 x is the coherent density of pumped polaritons in 
units of A x 2 . The diagonal process {k, k p } — > {k, k p } 
produces the renormalization of the energy dispersion 
(blueshift), namely E LP (k) = E LP (k) + 2V££ jV kp | 2 . 
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Notice the equality of the diagonal elements of M£ ar is 
equivalent to the exact energy conservation for the pro- 
cess {kp,k p } — > {k, kidier}- Finally, the coupling to 
the external photons is responsible for a radiative damp- 
ing 7^ = 2-Kg 2 \Ck\ 2 /h~ and a Langcvin force P k (i) = 
Jdttg C*e-* m a k ,n(0). 

From the coupled equations for the operators Pk(t) 
and f> k (t), we can obtain the physical quantities of in- 
terest. Let us consider the polariton occupation number 
Nk(t) = (Pk(t)pk(t)) ■ From the equations for p k (i), we 
get: 
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lm{(pl(t)F k (t)) 



In order to obtain the frequency spectrum of the para- 
metric luminescence, we have to come back to the equa- 
tions for pk(t) andp k (t). In the frequency space (through 
the transformation / °° e l ^ +l0+S)t dt), the coupled equa- 
tions read 



^k idler 



Pk(w) 

(2u p - 



M par 



Pk(w) 

(2u) v — LO) 



F k (w) + ihp^t = 0) 



i rpint 



Pi? e 



**<p£(*)pL.(*)>} 



This equation shows that the polariton population 
iVk(i) is driven by the anomalous quantum correlation 
(t)). The anomalous correlation evolution 



reads: 



i» S 04(*)pL,„(*)> = 
-^p(k)- J E iP (k idler )-i( 7fe + 7 feidler )] ( Pk (t)p kidler (t)) 



+ JVk(t) + \u„(t))) 



Notice that the emitter-idler correlation (p k (£)p kldler (£)) 
has a generation term proportional to the Boson enhance- 
ment factor (1 + -/V k (i) + ^V kldler (f)). This means that 
there is a spontaneous inhomogeneous contribution and 
a stimulation term due to the emitter and idler popula- 
tions. 



Neglecting the noise source F k (w) due to the external 
photons, the solution of this linear inhomogeneous sys- 
tem is 



A(w) ihpk(t = 0) 



- 2hujp. The 



with A(u>) = hui + £ L p(k id ie r ) + Hk, 
complex pole energies -E + k and -E_ k are the eigenval 
ues of the parametric matrix M^ ar and describe the 
new collective excitations of the microcavity system. 
Through the expectation value 7£e{(p k (w)p k (0))} = 
7£e{(p k (0)p k (c<j))}, we obtain the stationary lumines- 
cence as a function of the steady state population iV k = 
(pj c (0)p k (0)) and parametric correlation amplitude *4 k * = 
(Pk(0)pL,„(0)). The result is 
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FIG. 2. Contour plot of \E LP (k) + E L p(ki d i er ) -2E LP (k p )\ 
(meV) as a function of k = ^(sin# X , Sill Uy ).The 
white dashed line represents the zero value contour, 
i.e. exact energy- momentum conservation for the process 
{kp,k p } — ► {k, kidier}- The x-direction is that of the pump 
wave- vector k D . 



Simple algebra shows that steady-state solutions of the 
coupled equations for population and parametric correla- 
tion are ( Pk (t)p kidler (t)) = A s k e^ 2 ^* and JV k (t) = JV£. 
The anomalous correlation amplitude reads 



E int V 2 



IX. |2 _ (Tfc+^eJ 2 



\Ei nt Vi 



with<5 k = (2h\u p - E LP (k) ~ E L p(k idler )-i(-f k +"f kidle 
Moreover, the steady-state population is 
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^lm(Ei nt VlAtt). 
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Of course, -^ k . der = ^\Tlk /lki d ur- As we will show below, 
the PL signal ("as well as the parametric correlation am- 
plitude A^ and population iV k ) diverges when the pump 
density |P kp | 2 reaches a threshold value. In fact, above 
threshold, the pump polarization can be no longer treated 
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as a parameter. This means that the equation of motion 
for the pumped mode has to be included, accounting for 
the pump depletion. Moreover, when the transfer from 
the pump wave-vector to the other ones becomes very 
massive, multiple scattering can be no longer neglected. 
Hence, the results here presented are valid only below the 
threshold of the parametric luminescence. 
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FIG. 3. (a)Contour plot of the parametric photolumines- 
cence (log scale) as a function of the angle 6 X (deg) and energy 
(meV) for the pumped density n p — 0.03 n sat and 6 y — 0. 
The dashed-line is the unperturbed dispersion of the lower po- 
lariton branch, (b) Two-dimensional pattern of the spectrally 
integrated emission. 

For sake of illustration, we present our analytical re- 
sults by using GaAs parameters for a realistic microcav- 
ity structure with Rabi splitting equal to 7 meV. We 
consider the situation of zero exciton-photon detuning, 
i.e. E exc (0) — E cav (0). For simplicity, we consider a 
k-independcnt polariton linewidth 7^ w 0.4 meV (for 
large k, the decreasing of the radiative width is usu- 
ally compensated by nonradiative broadening). The two- 
dimensional dispersion of the lower polariton branch is 
shown in Fig. 1. We consider a pump which resonantly 
excites the polariton branch at the critical wave-vector 
such as Elp(0) + Epp(2k p ) = 2Epp(k p ), allowing the 
particular polariton fission {k p ,k p } — ► {0,2k p } (with 
increasing pumping, one has to consider the renormal- 
ized dispersion E L p instead of Elp). Fig. 2 contains the 
contour plot of the quantity \Elp(\s) + -Elp (kidier) — 
2E L p(k p )\ in k space (k p is along the x-direction) . The 
white-dashed line represents the zero value contour, i.e. 
exact energy conservation for the fission {k p ,k p } — > 
{k, kjdier}- Remarkably, the conservation is only weakly 
forbidden in a large squeezed portion of k-space (dark 



region) around k p . This implies that the parametric 
coupling is efficient on a broad angular range. This is 
actually found in Fig. 3 (a), where a typical contour plot 
of the photoluminescence (log scale) is shown as a func- 
tion of the angle 9 X (9 y = 0) and of the emission energy. 
Moreover, Fig. 3(b) shows the two-dimensional pattern 
of the spectrally integrated emission (linear scale). The 
depression of the emission at large angles is due to the 
vanishing cavity fraction |Ck| 2 . 
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FIG. 4. (a) Normalized emission lineshape around k=0 
(fix = 9y = 0; angular acceptance of 2 deg) for 5 densi- 
ties (A-E). Inset: log-log plot of the emission intensity versus 
pumped density (n sa t units), (b) Dispersions (meV) versus 
6 X (deg) for 6 y —Q at the highest pump density. Thin dashed 
line: unperturbed polariton dispersion Elp- Thick dashed 
line: blueshifted dispersion Elp- Thick solid line: lZe(E±) 
with IZe(E-) < TZe(E + ). The bifurcation occurs when the vi- 
olation of the energy conservation is large enough and implies 
negligible parametric emission. 

To determine the dependence on the pump intensity, 
we just need to analyze the evolution of the eigenval- 
ues -E±,k- For a given k satisfying exactly the energy- 
momentum conservation for the parametric scattering, 
the expression for the pole energies E±^ is very sim- 
ple, i.e. Tle{E ±M ) = E LP (k) and \Im{E± M )\ = 
(VlfeTfc.die. ± ^r'l^kpl 2 )- This means that the imagi- 
nary part of the pole energy -E-,k gets smaller at the ex- 
pense of the other pole E + ^. Consequently, the photolu- 
minescence spectrum PL(k, u>) gets spectrally narrower. 
The threshold density for the parametric luminescence is 
defined by im(£_ k ) = 0, i.e. E lnt \Vw p \ 2 = ^Jlklk idler - 
For this value, the PL signal, the parametric correla- 



tion amplitude and steady-state population N£ 



di- 



verge and therefore the pump depletion has to be in- 
cluded. All these features can be seen explicitly in Fig. 4. 
With increasing pump density the lineshape of the emis- 
sion around k — blueshifts and narrows (Fig. 4(a)). 
The output intensity (see the log-log plot inset) shows a 
threshold behavior as a function of the coherent density 
of pumped polaritons n p = |"Pk p | 2 /^x- Notice that when 
the polariton damping 7^ is small enough, the threshold 
density n t hr can be orders of magnitude smaller than the 
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exciton saturation density n sa t- With the realistic mate- 
rial parameters employed here, we have nthr ~ 0-05n sat . 
Finally Fig. 4(b), shows explicitly the energy disper- 
sions of the collective excitations. When the energy 
conservation for the parametric fission is well satisfied, 
TZe{E±{k)) ~ ElpQc). On the other hand, when the en- 
ergy conservation is strongly broken and consequently the 
parametric luminescence is negligible, TZe(E—(k)) and 
lZe(E + (k)) bifurcate, tending to the diagonal elements 
of the matrix M^ ar , i.e. £xp(k) and 2%ui p — -Exp(kidier)- 
In the intermediate case, consistently with experiments, 
the dispersion 7^e(£'±(k)) shows peculiar features such a 
change of sign of the in-plane group velocity |8| around 
2k p and a flattening |3]|| around k = (see also the 
emission spectra in Fig. 3). 

In conclusion, we have theoretically described the para- 
metric luminescence of microcavity polaritons, showing 
the subtle interplay between parametric correlation and 
stimulated scattering. The new collective excitations re- 
sulting from the anomalous emitter-idler coupling pro- 
duce a very peculiar luminescence pattern in frequency 
and momentum space, giving a key to recent experiments 
HHH- Further investigations are in progress for the 
regime above threshold not treated in this paper, but 



investigated in new experiments [|10|. We wish to thank 
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